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Abstract 

A consistent and intuitive description of the twist-4 corrections to the hadron 
structure functions is presented in a QCD-improved parton model using time- 
ordered perturbative theory, where the collinear singularities are naturally elimi- 
nated. We identify the special propagators with the backward propagators of par- 
tons in time order. 
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1 Introduction 

The twist-4 or inverse-power l/Q 2 corrections to scaling in deep inelastic scattering 
(DIS) processes are an important subject for a precision test of QCD and the analyses of 
deep inelastic scattering (DIS) data [1-4]. Jaffe and Soldate (JS) used OPE method to 
relate the inverse-power corrections to the target matrix elements of local operators [2] . 
On the other hand, the QCD-improved parton model can provide a simple and intuitive 
description of scaling phenomena as a compensate method for the OPE technique. In 
this aspect, Jaffe [4] recasted the results of work [2] in parton model. Ellis, Furmanski, 
Petronzio [3] (EFP) and Qiu [4] try to construct a complete QCD-improved parton model 
including the inverse-power short-distance corrections based on the manipulation of stand 
Feynman diagrams. Where EFP's diagram-approach needs a special complicate proce- 
dure to pick up the short-distance contributions. Qiu advanced the above technique and 
proposed that a set of special propagators can isolate the short- distance effect in Feynman 
diagrams [4]. However, the dynamic origin of the special propagators and collinear safe- 
ness of the theory are not yet clear. Obviously, we still lack a consistent QCD-improved 
parton model to analyze the inverse-power corrections in DIS processes. 

In this work, we try to realize the twist-4 short-distance corrections to the hadron 
structure functions using time ordered perturbative theory (TOPT). We find that the 
collinear singularity can be naturally eliminated if we generalize TOPT to the collinear 
limit. In this framework, the twist-4 short-distance corrections to the hadron structure 
functions will present an intuitive space-time picture. In particular, we identify Qiu's 
special propagators [4] with the backward propagators of partons in time order. Our new 
method provides a consistent way to study higher-twist processes in a QCD-improved 
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parton model. 

The outline of the paper is as follows. In Sec. 2, we shall generalize TOPT to 
the collinear limit. We find that the collinear and infrared singularities in amplitude 
will transfer to the interaction-time from the momentum-space in TOPT and they can 
be natural eliminated. Using this TOPT-framework, we shall give a set of calculating 
rules, which pick up the inverse-power short-distance corrections to the hadron structure 
functions in Sec. 3. We shall discuss the factorization scheme in Sec. 4. A summary is 
given in the last section. 

2 Collinear safeness 

As we know, TOPT is a natural way for studying space-time properties of the twist-2 
DGLAP evolution equation [5]. In principle, a Feynman diagram can be decomposed to 
all possible time-ordered graphs in TOPT. However, in the leading-twist processes many 
TOPT amplitudes, which contain the backward-propagating partons, are suppressed in 
the infinite momentum frame (see Fig. 1). The reason is straightforward. The cross 
section of DIS is generally a function of the longitudinal momenta, x^P and transverse 
momenta, l\, of partons. The terms with x { P are canceled in Fig. la due to the con- 
servation of longitudinal momenta at every vertex and the final results depend only on 
Ij 2 and give log Q 2 -corrections. Any backward-propagating parton in the splitting vertex 
will violate the balance among XiP in the energy terms and lead to dl\j P 2 — > 0. Thus, 
the backward-propagating partons are suppressed and the scattering amplitudes only con- 
tain the forward-propagating partons in the leading twist processes. In this work we find 
that the twist-4 amplitudes also expose the simple TOPT-form but the contributions of 
backward moving partons are no longer neglected. 
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The higher twist corrections to structure functions are separated in the collinear ex- 
pansion, where all transverse momenta of partons are taken as zero [4]. Therefore, we 
first generalize the TOPT-amplitude to the collinear limit. As an example, we consider 
the process p\ + p 2 — > p[ + p' 2 , in which the parton momenta are parametrised as 



Pi ■ (PW,l±,Plz), P2 ■ (P70,-1±,P2z), P'l ■ (PlO,l±,Plz), P 2 : (PlO > ~ l ±, P2 z) > k '■ (k ,0±,k z ) 

(1) 

where the transverse momenta l± will be set to zero. We emphasize that our following 
results are independent of the parameters (1). We choose a coordinate frame such that 



p n = pn^, 

q» = -x B plf + ■^—n". (2) 
2x B p 

The contributions of the Feynman propagator of fermion (quark or antiquark) with 
the momentum k to the amplitude are 



= — { / dt i e -it'{u-p' w -p' 20 ) / dt e -«(Pio+P20-o;)[_ i7 . k6 ( k(j _ + 
2cJ J —oo J —oo 

r dte- lt{pw+P2n+u) f ^ / e i *'Wo+f2o+^)[i 7 . fc«y(jfe + to)]} 

J—oo J —oo 

(2n) 3 5 3 (p 1 +p 2 -0 1 -0 2 ) 
= [S F (k) + S B (k)](27r)W(p 1 +p 2 -p' 1 -p 2 ), (3) 
where to = \Jk\ + (k z ) 2 . We use light-cone coordinates, in which we define 
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n^ = (n°,n ± ,n 3 ) = ^=(l,0 ± ,l); n" = -^=(1, ± , -1), (4) 



or equivalently 



rF=(ri+,rr,n ± ) = {l,0,0±y, n» = (0, 1, ± ). (5) 

Thus, n 2 = 0,n 2 = 0, and n • n — 1; x — k • n/P ■ n is the light-cone fraction. In the 
collinear limit, we take l± — > 0. In consequence, we have 



]imS F (k)(27r)W(p 1+ p 2 -ft-p' 2 ) 
= lim 7-M(fe -q;) r~ dt>e - it ^- P ' 10 - P ' 2o) 

p-it' (pio +P20 -u>) (3 -«7 1 (pio+P2o - w) 

(^— )(2vr) 3 5 3 (p 1 + p 2 - £ - 

PlO + P20 - W pio + P20 - W 



2cu 



- / dt'e^'^-^-^it 1 - T)(27r) 3 5 3 (p 1 + p 2 - ft x - ft 2 ) 

J — oo 



\f2a^ ■ n 



(2n) 4 6\p 1+ p 2 -ft-p' 2 ). (6) 



2cj 

Usually, the lower limit T of the integral, which is relating to the interaction time, is 
taken as — oo. However, T is really a finite quantity in physics and we denote t'—T = a/u, 
where a is a dimensionless quantity. We emphasis that the term — e _lT ^ 10+P20_ ^/(Pio + 
P20 —oj) in Eq. (6), which vanishes in the noncollinear case, is important for the collinear 
safeness. 

The second term in Eq. (3) is 



lim . S B (k)(2n) 4 5\ Pl +p 2 -ft-p> 2 ) 

£j_ — ►U 
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P'lO+P'20^ 2CJ J -00 Pl0+P20+ W - - 

( ) 5 (Pi +P2 "Pi -P 2 ) 



4w 2 

= ^(2 7 r) 4 5 4 ( Pl+P2 -K- P ' 2 ). (7) 

The negative frequency solution fc = —u in Ssik) can be understand as a backward- 
moving parton in time order since k5(ko+u) = —(a;, 0j_, — & 2 ) = fc, fc is on-shell momentum 
in TOPT. We will refer to Sp(k) and Ssik) as the forward and backward propagators, 
respectively. We can generalize k to the off-shell momentum k in (6-7) since n ■ n — 0, 
where the transverse momentum of parton is vanished. Obviously, the backward propaga- 
tor of quark is identical with the special propagator [4] or the instantaneous part of quark 
propagator in the non-collinear case [6] . Qiu has pointed out that the special propagator 
offers no space separation along the light-cone coordinate of two interaction points con- 
nected by the propagator. One can straightforwardly understand the above mentioned 
property as follows: the propagation of a backward parton "seems" to be instantaneous 
in kinematics when P — > 00. It is interesting that Fig. 2b shows that the virtual photon 
creates a qq pair in the infinite-momentum frame. 

In general, the massless partons with the parallel momenta can go on-mass-shell simul- 
taneously in the collinear case and collinear singularities may arise. Therefore, one might 
think that divergence could arise in the forward propagator Sp(k) if k 2 = (we call it as 
the collinear propagator). Fortunately, Eq. (6) shows that the collinear singularity is ab- 
sent due to the compensative contributions from the factor — e - jT (pi()+P2o-w) j +p 20 — cj) 
and a/uj is a finite quantity in Eq. (6). Furthermore, due to helicity conservation, the 
contributions of the numerator containing the collinear forward propagator to the twist-4 
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structure functions vanish. For example, the emission of an on-shell massless collinear 
gluon along the quark line is forbidden in the collinear approximation. Therefore, one 
can neglect the collinear propagator Sp(k) in the twist-4 component of structure function 
provided a/co is finite. Thus, we show a new way to eliminate collinear divergence: the 
collinear singularities will be transferred from the momentum-space to the interaction- 
time in TOPT and the later is really finite in physics. 

Now let us discuss the gluon propagator G af3 (k). Similar to Eqs. (3-7), we have 



lim / d 4 kG af3 (k) 



2aT^(k)5(k -u) ^ 4 4 , T^(k)5(k + u) ^ 4 4 , 
= (kTn)~ 2 ( ^ (P1+P2-P1-P2) ^r—j2 (2tt) 5 (P1+P2-P1-P2) 

ee [Gf{k) + Gf{k)]{2nY5\ Pl +p 2 - pi -p' 2 ), (8) 
where T alS is the gauge-dependent polarization sum. In the light-cone gauge and limit 



T^(k)5(k -u) = -g aP + ^ + kf>na 6(ko - uj) 

Zco 



= (9) 

in Gf{k) 
and 

h a nP -I- hPn a 

-T^(k)6(k + u,) = <?* - n o +fc U 5(k + co) 

2u) 

in G^(k), respectively. Where we use to = xp. Therefore, 



and 



where 5^ and g a/3 collect the contributions of the terms with l\ in the numerator. Thus, 
in the collinear limit, we need only consider the contributions of 



Gf(k) = (13) 



i a n 

in calculations of the twist-4 structure functions. Equation (13) is the special propagator 
G"P(k) [4] or the instantaneous propagator for gluon in the noncollinear case [6]. 

From the above discussions, we can conclude that the contributions of the vertices 
(Figs. 3a and 3b) vanish due to collinear safeness, where F and B label the forward- and 
backward-moving partons, respectively. Furthermore, the vertices in Figs. 3c and 3d also 
are forbidden due to 7 • ny ■ n = and n • n — 0. 

The amplitude with twist-4 corrections in the collinear expansion can contain the 
noncollinear propagators with off-shell momenta although parton transverse momenta are 
neglected. Any noncollinear propagators can be decomposed to a non-vanishing forward 
and backward components in TOPT, in which Ss{k) and G% (k) have the same form as 
(7) and (13), respectively. On the other hand, we can use stand TOPT to calculate S F (k) 
and Gp{k) and get 



(14) 



and 



Gf(k) = j-. — — : r, (15) 

where k are off-shell momenta. From the discussions about the propagators, we can 
find that the backward-propagating gluon is longitudinally polarized, while the forward- 
propagating gluon is transversely polarized in the collinear approximation. A parton will 
inverse (or will keep) its moving direction after absolving a forward (or backward) gluon 
due to helicity conservation (Fig. 4). It means that the vertices of Figs. 3e-g are forbidden. 
Thus, a backward-moving gluon only couples with two forward-moving partons, while a 
forward-moving gluon couples with a pair of quark and antiquark moving oppositely. 

The selecting rules (Fig. 4) simplify the calculations of the twist-4 corrections in the 
collinear expansion. In fact, we have the following TOPT-procedure for calculating the 
high-twist contributions in the collinear expansion: (1) to list all possible time-ordered 
diagrams of a photon-parton forward scattering amplitudes; (2) to label B or F on the 
parton lines according to Fig. 4, and note that all loop momenta are forward while 
the cut propagator is denoted by F or B for gluon or fermion, respectively, since the 
polarization sums of initial and final gluons have the same form S^; (3) the contributions 
of the backward and noncollinear forward propagators are calculated by Ssik), G^(k), 
S F (k) and Gp^(k), respectively. One can find that every cut-diagram only corresponds 
to a non- vanished TOPT-graph due to the coupling of partons in space-time (Fig. 4). 
In opposition to our TOPT-method, the covariant perturbative methods [3,4] do not tell 
whether forward propagator or backward propagator is more important, since they mix 
together all possible time ordering in Feynman diagram. Obviously, the gauge invariance 
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is preserved if we include all TOPT-diagrams of a physical process at given order of 
coupling constant. 

3 Power corrections to structure functions 

One can simply reproduce all results of Ref. [4] using our TOPT-programs in Sec. 2. 
In deep inelastic scattering of lepton-hadron, the hadronic tensor W^ v (p, q) can decompose 
to 

W^(p, q) = e?(x B , Q 2 ) + eF(x B , Q 2 ), (16) 

where and are the longitudinal and transverse tensors, respectively. Using the 
definition (4), 



I (, r + l£g>W + - 1(*V)) , (17) 



= 1 (n*W + n^rf - g^) = \d? v . (18) 
The standard structure functions are defined as 

2F 1 (x B ,Q 2 ) = F T (x B ,Q 2 ), 

F 2 (x B , Q 2 ) = x B [F T (x B , Q 2 ) + F L (x B , Q 2 )}. (19) 

The inverse-power corrections of the twist-4 QCD processes to the structure functions 
F\ and F 2 , or the hadronic tensor are from Fig. 5 (for quark-gluon-to-quark-gluon) 
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and Fig. 6 (for two-quark-to-two-quark), respectively. In the concrete, W£ v ~ / a^T, 
where a 1 *" are the cross section of the processes in Figs. 5 and 6, T is the contributions 
from the target part [4,7]. Using our calculating program in TOPT, we can easy write 
<t^ u in Fig. 5 as 



^qga( X , X U X 2 ) 

= 2n^5(x-x B )(^-) 2 , (20) 



(j%g b (X,X 1 ,X 2 ) 



167T 



Tr[Y 



X BP _ a Q 2 R X BP _ „ _ 1o , ( ,XB 

= 2n^5(x-x B )(^) 2 , 



(21) 



d a/3 Q 2 a x B p _ 7-n v . xb 

i6^ Tr[y W • n7 q^ 7 • " Afr-^sp 7 P7 ' n]27r5(a;2 " Xb) -q 2 



= e^8(x 2 - x B ) 



x 2 



(xi - x 2 )Q 



2 ' 



(22) 



Vqqe{ X , X li X -2) 



16tt 



Tr[j' 



, x BV^ Q 



2 



q* -7 ' n7 "2^ 7 ' nr 2*7 7 ^ 7 • n]27r5(x - Xb)X 4 

= rtn v b{x - x B )^-. 

Q 



(23) 
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Note that the contributions of Figs.5g-5j are vanished since in Figs. 5g and 5h the 
propagators Fi contributes ~ 7 • n/p and leads to o^ 1 jp 2 — > if p — > 00; Figs. 5i and 5j 
include Tr [7^7 • ny a y ■ ny v y ■ ny^j • n] — 0. 

Similar, the contributions of four-quark process to the leading inverse power corrections 
in Fig. 6 are 



16tt 



7 p7 • nj 



2(x 2 - x B )p 

= ej^8{x — x B 



— —7 p7 • n \27r5(x - x B ) — 

2{x 1 -x B )p Q 2 



x 



(x 1 - x)(x 2 - x)Q 2 ' 



(24) 



^ rr r « 7 - ™ a - /3 a -i 2;r l or/ \ X ^ 

lfa rrh 2(x 2 -x B )p 7 " ' m 2( X2 - x B )i ' " ' n] T^W 2 * Hx > ~ XB) Q 



2 



= ^ S(xi - XB) j —^—^ i . (25) 
Thus, we reproduce the results of Refs. [3,4] using TOPT. They are 



a%(x,x 1 ,x 2 ) 



JL(e!? + e£)6(x - x B ) 



2x B 



5(x 2 - x B ) - S(xi - x B ) 



x 2 - x 1 



(26) 



, S(y 2 



5(x 2 - x B ) — 8(x 1 — x B ) 5(x - x B ) - 5(x 1 - x B ) 



(x 2 - x)(x 2 - Xi) 
xb)-S(vi-xb) S(y 



(x 2 — x)(x — Xi) 
x b) -S(yi- x B 



(V2 ~x)(y 2 -yi) 



{y 2 -x){x-y 1 ) 



(27) 
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where y\ = x\ — x and 1/2 = x — x 2 . 

In the above mentioned calculations we considered the processes of two-parton-to-two- 
parton. We also should pick up the contributions from all cut diagrams of the interference 
process of one-parton-to-one-parton and three-parton-to-one-parton. For example, we 
need include contributions of Fig. 7, which are corresponding to Fig. 5c. One can find 
that the contributions of Figs. 7a and 7b are the same as that of Fig. 5c. Therefore, we 
should multiply the above mentioned results (26) and (27) by a factor 3 for including the 
contributions of the interference processes. 

4 Factorization scheme in TOPT 

A relating subject with the above mentioned statements of power corrections to struc- 
ture functions is the factorization scheme in TOPT. As we know that factorization theorem 
allows us to separate the calculable short-distance effects from incalculable long-distance 
part. A complete derivation of the factorization theorem in DIS has been given by many 
authors [8]. The discussion of this section intends to provide an intuitive understanding 
of the existing factorization scheme. 

A simplest amplitude of DIS process contains a pair of (left and right) loop Feynman 
propagators and it can be written as 

M^ip, q) = Tr[ 7 • pTDHNH*D*T*] , (28) 

where D and D* are Feynman propagators of quarks connecting with the target-part T 
(T*) and hard-part H (H*), respectively. According to the view of parton model, in a 
factorization scheme, all virtual (off-shell) partonic lines joining T and H can be broken 
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to the free (on-shell) partonic lines. In this case, we can re-stipulate those on-shell partons 
as incoming and outgoing partons in the target- and hard-parts, respectively. 

Any Feynman propagator is decomposed to the forward F and backward B propaga- 
tors and the propagating partons are evaluated on-shell in TOPT. The diagram including 
two forward propagators in Fig. 8a gives the leading twist-2 contribution and we can 
separate the calculable short- distance effects from incalculable target part, since two on- 
shell-forward propagators can be broken to incoming and outgoing partons in the target 
and hard parts in a factorization scheme, i.e., 

M^(p, q) = Tr[ 7 • pT 7 • kHNH*^ ■ kT*} 

= Tr[T* 7 • pT 7 • k]Tr[HNH*-f ■ k], (29) 

where k is on-shell momentum of the propagator F. The first factor in Eq. (29) is the 
trace of the (soft) forward virtual quark-target scattering amplitude, and the second factor 
relates to the (hard) probe-parton subprocess. 

On the other hand, the backward propagator implies the contact (or instantaneous) 
correlation and we can not re-stipulate the backward propagator in any factorization 
scheme. Therefore, one of the factorization conditions is that we can avoid the backward 
propagators between T and H. For this end, we note that a backward parton always 
couples with a pair of forward partons according to Fig. 4. After "pulling down" two 
forward propagators from the target matrix, the backward propagator can be isolated into 
a hard-subgraph of a next twist amplitude. In consequence, we get the twist-3 (Figs. 8b 
and 8c) and twist-4 contributions (Fig. 8d), respectively. In this case, all initial partons 
are forward moving along the direction of target and we can break and recombine those 
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parton lines between two factorized parts in the factorization scheme. Through the above 
TOPT-expansion, we also show a clear picture where the different twist contributions are 
separated. 

5 Summary 

In this work, we presented following physical picture about the inverse-power correc- 
tions to the hadron structure in parton model: (i) a hadron is assumed to be consisted 
by forward initial partons along the target momentum in the infinite momentum frame 
(i.e., the positively of initial parton momentum) and the transverse momentum of parton 
can be neglected if comparing with its longitudinal momentum; (ii) a probe strikes with 
a parton and forms a covariant Feynman propagator connecting with the target matrix 
and probe; (iii) the Feynaman propagator decomposes to the forward and backward prop- 
agators in TOPT and the backward propagator "pulls down" two forward propagators 
from the target matrix; (iv) after factorization, two-parton-correlation function has the 
dimension of [length]' 2 , therefore, the hard-amplitude of two-parton-to-two-parton will 
provide a contribution with 1/Q 2 to the total amplitude. This leads an inverse power 
correction to the hadron structure function. 

In summary, we have generalized time-ordered perturbative theory to the collinear 
limit, where the collinear singularities are eliminated in a natural way. We find that 
a clear space-time picture of the twist-4 amplitude is presented in its TOPT-form. In 
particular, we identify the special propagators with the backward propagators of partons. 
Our method provides a new way to analyze the twist-4 processes in a consistent QCD- 
improved parton model. 
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Figure Captions 

Fig. 1 A twist-2 photon-hadron scattering process of 0(a s ) in TOPT; (a) dominant 
amplitude with the forward propagator and (b) suppressed amplitude with the backward 
propagator. 

Fig. 2 A twist-4 photon-hadron scattering process of 0(a s ) in TOPT; (a) suppressed 
amplitude with the forward propagator and (b) dominant amplitude with the backward 
propagator. 

Fig. 3 Some forbidden vertices of partons in twist-4 photon-hadron scattering amplitude. 
F and B indicate forward- and backward-moving partons, respectively. 
Fig. 4 Elemental vertices in the collinear approximation. 

Fig. 5 A complete set of quark-gluon-to-quark-gluon TOPT diagrams, which give leading 
power corrections to the hadron structure function. 

Fig. 6 Non-vanished two-quark-to-two-quark TOPT diagrams, which give leading power 

corrections to the hadron structure function. 

Fig. 7 Interference processes corresponding to Fig. 5c. 

Fig. 8 TOPT-expansion of the different twist contributions. 
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